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ABSTRACT 

Observations show that early-type galaxies contain a considerable amount of interstel- 
lar dust, most of which is believed to exist as a diffusely distributed component. We 
construct a four-parameter elliptical galaxy model in order to investigate the effects of 
such a smooth absorbing component on the projection of kinematic quantities, such 
as the line profiles and their moments. We investigate the dependence on the optical 
depth and on the dust geometry. Our calculations show that both the amplitude and 
the morphology of these quantities can be significantly affected. Dust effects should 
therefore be taken in consideration when interpreting photometric and kinematic prop- 
erties, and correlations that utilize these quantities. 

Key words: dust, extinction - galaxies : elliptical and lenticular, cD - galaxies : ISM 
- galaxies : kinematics and dynamics 



1 INTRODUCTION 

Theoretical work on the transfer of radiation in the inter- 
stellar medium demonstrated that the effects of dust grains 
are dominant over those of any other component (Mathis 
1970; Witt et al. 1992, hereafter WTC92). From the ob- 
servational point of view, dust in elliptical galaxies is opti- 
cally detectable by its obscuration effects on the light dis- 
tribution, when it is present in the form of dust lanes and 
patches (e.g. Bertola & Galetta 1978, Hawarden et al. 1981, 
Ebneter & Balick 1985, Sadler & Gerhard 1985). In emis- 
sion, cold dust [Td ~ 30 K) is detected at far-infrared wave- 
lengths. The /i? AS satellite detected dust in ellipticals in the 
60 and 100 /xm wavebands (e.g. Jura 1986, Bally & Thron- 
son 1989, Knapp et al. 1989) in, at that time, unexpectedly 
large quantities. However, IRAS is insensitive to very cold 
dust {Td < 30K), and preliminary ISO results (Haas 1998) 
and (sub) millimeter data (Fich & Hodge 1993, Wiklind & 
Henkel 1995) may indicate the presence of an additional very 
cold dust component in ellipticals. 

All these observations show that the presence of dust in 
ellipticals is the rule rather than the exception (van Dokkum 
& Franx 1995, Roberts et al. 1991). A straightforward ques- 
tion is whether the dust detected in emission and absorp- 
tion represents one single component. Goudfrooij & de Jong 
(1995) show that the IRAS Aust mass estimates are roughly 
an order of magnitude higher than those calculated from 
optical data. To solve this discrepancy, they postulated a 



* Research Assistent of the Fund for Scientific Research 
ders (Belgium) 



Flan- 



two-component interstellar dust medium. The smaller com- 
ponent is an optically visible component in the form of dust 
lanes, the larger one is diffusely distributed over the galaxy, 
and therefore hard to detect optically. At least part of this 
smooth dust distribution can be easily explained in ellip- 
ticals : the old stellar populations have a significant mass 
loss rate of gas and dust (Knapp, Gunn & Wynn- Williams 
1992), which one expects to be distributed smoothly over 
the galaxy. 

The effects of a diffuse dust component on the photom- 
etry are studied by WTC92 and Wise & Silva (1996, here- 
after WS96). They demonstrate that even modest amounts 
of dust can produce significant broadband colour gradients, 
without changing the global colours. Radial colour gradi- 
ents can therefore indicate the presence of this diffuse dust 
component. 

Hitherto the effect of a diffuse absorbing medium on the 
internal kinematics of ellipticals has been neglected, proba- 
bly due to the widespread idea that early-type galaxies are 
basically dust-free. In this paper, we construct a set of el- 
liptical galaxy models in order to investigate these effects. 
Our models contain 4 parameters, allowing us to cover a 
wide range in orbital structure, optical depth and relative 
geometry of the dust and the stellar component. 

We want to demonstrate the qualitative and quantita- 
tive effects of dust absorption on the observed kinematic 
quantities. In the stellar dynamical view of a galaxy, the ul- 
timate purpose is the determination of the phase space dis- 
tribution function F(r,v), describing the entire kinematic 
structure. This distribution function is usually constructed 
by fits to observed kinematic profiles as the projected den- 
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sity pp{x) and the projected velocity dispersion ap{x). Us- 
ing a completely analytical spherical model for ellipticals, 
we demonstrate in this paper the effects of dust on these 
projected profiles and on line profiles. 

In section 2 we describe the mechanism of radiative 
transfer and the effect of dust on the projection of physi- 
cal quantities. In section 3 we describe our model in detail. 
The results of our calculations are presented and discussed 
in section 4. Section 5 sums up. 



2 DUSTY PROJECTIONS 

A fundamental problem in the study of stellar systems 
is to retrieve three-dimensional information from a two- 
dimensional image. This deprojection problem is degenerate 
and can only be solved with specific assumptions, e.g. spher- 
ical symmetry. In this section we show how dust absorption 
affects the projection. In this first study we limit ourselves 
to absorption effects and we do not include scattering. 

The transfer equation describes how radiation and mat- 
ter interact. The radiation field at a point r propagating 
in a direction n is characterized by the specific intensity 
I{r, n), while the global optical properties of the galaxy are 
described by the light density p{r) and the opacity K(r). 
The transfer equation on a line-of-sight p is then 



observer 



■j^(p;s) = p{r) 



K,{r)I{p;s) 



(1) 



where s is the path length on p and I{p; s) = I{r, Ss) is the 
intensity at the point s on p propagating in the direction of 
the observer e^. Equation (hi) can be readily solved to obtain 
the light profile detected by the observer pp{p) = I{p; So), 



Ppip) 



p{r) exp 



K(r) As' ) ds 



(2) 



where So denotes the location of the observer and the inte- 
gration covers the entire path. Equation (g) can easily be 
generalized to include kinematic information: 



fp{P,v) = / .f{r,v) exp 



K(r) ds' I ds 



(3) 



This expression yields what we call the dusty projection 
fp{p, v) of the kinematic quantity /(r, v). 

To describe the observed kinematics of spherical galax- 
ies we adopt a coordinate system as illustrated in figure hia. 
The distance between the observer and the centre of the sys- 
tem is denoted by D. Due to the spherical symmetry we can 
write f{r,v) — f{r,v) and «;(r) = «;(r), with r — \\r\\ the 
spherical radius. Lines-of-sight are determined by the an- 
gle t/j they make with the line-of-sight through the system's 
centre. Substituting this in equation (H), one finds 



fpi^,v) =exp 

f-D 



D 



K{r) r dr 



X 



exp 



'r^ — x^ 
K{r')r' dr'\ f{r,v)rdr 



V^ 



+ 



exp 



K(r')r'dr'\ f{r,v)rdr 



(4) 



\/r^ — x'^ 

where we have set x = D sintp. 

This equation can be simplified substantially by mak- 
ing the approximation of parallel projection, as in figure nib 




Z = +00 observer 



D 



z = 



(a) 



(b) 



Figure 1. The two coordinate systems as described in section 2. 
When the distance D between the observer and the galaxy's cen- 
tre goes to infinity, figure (a) transforms into figure (b). 



Since the scalelengths of galaxies is several orders of mag- 
nitude smaller than their distance, the induced error of the 
order (x/D)'^ is usually negligible. In this system a line-of- 
sight p is determined by its projected radius x on the plane 
of the sky. Equation (pi) becomes 

/-t-oo / P + CG \ 

f{r,v) expl- K{r)dz'jdz. (5) 

We can rewrite the integrals in (^j as integrals over r and 
obtain 



fp{x,v) 


J X 


rCr r)-"- ' ' 
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where 






IC{x,r) 


= exp f - 


/■+°° K{r)rdr\ 

X coshfr'^('->'d^' 



(6a) 



(6b) 



The same result can sirnply be found by taking the limit 
D — > -l-oo in expression (14). 

The weight function IC{x,r) assumes values between 
and 1, and reduces to unity if the opacity vanishes, in cor- 
respondence with the normal, dust-free projections. In the 
sequel of this paper we will always refer to parallel projec- 
tion. 



3 THE MODEL 

3.1 The stellar component 

Since dusty projections already require a double integration 
along the line-of-sight, one is inclined to use a dynamical 
model that is as simple as possible to avoid cumbersome 
numerical work. Therefore the Plummer model (Dejonghe 
1987), representing a spherical cluster without central sin- 
gularity, seems an obvious candidate. Although the photom- 
etry of this model does not fit real elliptical galaxiesn, it 

* A more thorough analysis of the effects of diffuse dust on more 
realistic photometric profiles can be found in WTC92 and WS96. 
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has the huge advantage that the entire kinematic structure 
is sufficiently rich to include radial and tangential models, 
and that it can be described analytically. We can safely as- 
sume that our results will be generic for the class of elliptical 
galaxies. 

We fix the units by choosing a value for the total stellar 
mass M, the gravitational constant G and a scale length 
r, . Working in dimensionless units we choose M = G = 
r« = 1. For the potential we choose the natural unit ?/), = 
GM/rt and we set V(0) — 1- In these units the central 
escape velocity equals \/2, and the Plummer model is defined 
by 



ij{r) 



VT+r^ 



^W = i^F (1 



(7a) 
(7b) 



Dejonghe (f986) describes a technique to construct 2- 
integral distribution functions F{E, L) that self-consistently 
generate a given potential-density pair, based on the con- 
struction of an augmented mass density p(^/', r), i.e. the den- 
sity expressed as a function of r and V'- For the potential- 
density pair (Falb) the following one-parameter family is 
convenient 



p{'^,r) = j^^^ 



1 + r 



(8) 



This results in a family of completely analytical anisotropic 
Plummer models, fully described by Dejonghe (1987). The 
distribution function Fq(E,L) and its moments can be ex- 
pressed analytically. For example, the radial and tangential 
velocity dispersions read 



alir) = 



1 



1 



(^^r) =(Jg{r) = 



1 



6-gyi~H 



1-i. 



2 1+r^ 



and Binney's anisotropy parameter equals 



/3(r) = 1 



2 l + r2 



(9a) 



(9b) 



(10) 



This expression provides the physical meaning of the pa- 
rameter q. It has the same sign as I3{r), and its range is 
restricted to — oo < g ^ 2 (since — oo < I3(r) ^ 1). Tangen- 
tial clusters have q < Q, radial clusters have < g ^ 2 and 
for q = Q the model is isotropic. For small radii all models 
are fairly isotropic, while the true nature of the orbital struc- 
ture becomes clearly visible in the kinematics at radii r 2> 1. 
More specifically, we consider a tangential model (g = —6), 
a radial one {q = 1) and the isotropic case [q — Q). 



3.2 The dust component 

A dust model is completely determined by its opacity func- 
tion fv(r). Besides the geometrical dependency, this function 
is also dependent upon the wavelength, which we do not 
mention explicitly. WTC92 considered an elliptical galaxy 
model with a diffusely distributed dust component. They 
used a King profile for both the stellar and the dust compo- 
nent. Analogous models were considered by Goudfrooij & de 
Jong (1995) and WS96. In the assumption that dust is only 
present in the inner part of ellipticals, WTC92 cut off the 



dust distribution at two thirds of the galaxy radius. Never- 
theless dust may also be present in the outer regions, with 
temperatures of the order of 20 K (Jura 1982, Goudfrooij & 
de Jong 1995). 

Our dust model is defined by the opacity function 



<r)^'r'^-^ 



^ r (^) ro 



lfl + ^ 



(11) 



with ro a scale factor and a the so-called dust exponent. The 
normalization is such that r equals the total optical depth 
(in the V-band), defined as 

/"-f oo 

r = 2 / K(r)dr. (12) 

Jo 

We can calculate the weight function K,{x, r) by substitut- 
ing (0) in (©, 



K,(x,r) — exp ( \—!- ) cosh 



t{x) 



U„- 



J.2 _^ J.2 



.(13) 



where r(x) is the optical depth along a line-of-sight x 



-{x) = 2 



K(r)rdr / x 



(14) 



and the function Uq (z) is defined in appendix ^. 



3.3 Parameter space 

3.3.1 The optical depth 

The optical depth is an indicator of the amount of dust 
in the galaxy, and is therefore the primary parameter in 
the models. WS96 estimated the optical depths of a sample 
of ellipticals by fitting radial colour gradients and found a 
median value of r ;$ Sri We explore optical depths between 
r = 0.1 and r = 200. 



3.3.2 The dust geometry 

The parameters ro and a determine the shape of the dust 
distribution. The exponent a is restricted to a > 1. Smaller 
values of a or larger values of ro correspond to extended 
distributions; for larger values of a or smaller values of rp 
the dust is more concentrated in the central regions. Since 
the two geometry parameters ro and a have roughly the 
same effect, we concentrate upon the case ro = 1 and we 
take a as the parameter determining the geometry. 

For high values of a the distribution is very concen- 
trated in the central regions of the galaxy. HST imagery 
has revealed compact dust features in the majority of the 
nearby ellipticals (van Dokkum & Franx 1995), and these 
could be an indication of an underlying concentrated diffuse 
dust distribution. However, Silva & Wise (1996) showed that 
such smooth dust concentrations should induce easily de- 
tectable colour gradients in the core, even for small optical 
depths. But high resolution HST imaging of two samples of 

* This is scaled to the definition of r adopted in this paper, de- 
fined as the projection of the opacity along the central line-of-sight 
through the entire galaxy. WS96 defined r as the integral of the 
opacity from the centre of the galaxy to the edge, half our value. 
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nearby ellipticals revealed that elliptical galaxy cores have 
usually small colour gradients (Crane ei al. 1993, CaroUo 
et al. 1997), and thus no direct indication of concentrated 
diffuse dust is found. Therefore high values of a seem to be 
less probable and we will only treat them for the sake of 
completeness. 

On the other hand, for small values of the dust com- 
ponent the distribution will be extended. If a approaches 1 
the distribution is very extended, and little dust is located 
in the central regions of the galaxy. In the limit the dust is 
infinitely thin distributed over all space, the contribution of 
the dust in the galaxy is zero and the dust effectively forms 
an obscuring layer between the galaxy and the observer. This 
geometry is known as the overlying screen model (Holmberg 
1975, Disney ei al. 1989) and is the analogue for extinction 
of galaxian stars. By taking the limit q ^ 1 in (n4) we 
obtain 



t{x) = T 

and substitution in (hsl) yields 



(15) 



/C(a;,r) =exp(^--j . (16) 

Notice that the same result can be found by taking the limit 
ro -^ +CO. 

As demonstrated in section 2, the computational prize 
of dust absorption on the projection of kinematic quanti- 
ties is a double integration along the line-of-sight instead 
of a single one. This can be avoided if the weight function 
lC{x,r) can be calculated analytically for all x and r, i.e. 
if the functions 11^-3(2) can be evaluated explicitly. For 
obvious computational reasons, we concentrate on the cases 
where a is an integer number. But to study the effects in the 
neighborhood of the singular value a = 1, we also consider 
a range of non-integer values for a between 1 and 2. 



^.1.1 The projected density 

The primary effect of dust absorption on the projected den- 
sity profile (figure pk,b) is obvious : an overall reduction. 
Even for modest amounts of dust this effect is clearly ob- 
servable. The projected density in the central regions is al- 
ready reduced to half of its original value for an optical depth 
r = 1.5. At the core radius Tc — 0.64 this occurs for r = 3. 
The global effect of the extinction is shown in figure Be, 
which displays the observed integrated light Mobs, i-e. the 
fraction of the total light output that can be detected by the 
observer. It is obtained by integrating the projected density 
profile Pp{x) over the entire plane of the sky. 



Afohs = 27r 



Pp{x) xdx. 



(19) 



The global extinction effect is strong : for an optical depth 
r = 1 already 20 per cent of the mass is hidden for the 
observer, and for r = 3.4 this becomes 50 per cent. 

Besides reducing the amount of starlight, the dust has 
also a qualitative effect on the projected density profile. 
Since the dust is primarily located in the central portions 
of the galaxy, the pp{x) profile will be flattened towards the 
centre. As shown by WTC92 and WS96 this causes the pro- 
duction of colour gradients. But it has also a non-negligible 
effect on all global photometric quantities, e.g. luminosities 
and core and global diameters. For example, the apparent 
core of the galaxy (quantified by the core radius r^ or the 
effective radius r^) will be larger due to the central flatten- 
ing (figure |2|d). This effect can be up to 25 per cent or more 
for our basic model with moderate optical depths t « 2. 
On the other hand the total diameter of the galaxy (usually 
expressed as the diameter of the 25magarcsec~^ isophote) 
will appear smaller due to the global attenuation. 



4 THE PROJECTED KINEMATICS 

For the dust-free model all the moments of a line profile can 
be calculated analytically (see appendix y). For the lower 
order moments we find 



(17) 



''^^'^ -^(i+x^r 






and 






2, ^ 9vr 1 1 

(T I T 1 — 


('- 


n 2 

bq X 


""^^""^ 32 6-gVl + x-2 


12 1+x^ 



(18) 



For dusty models this projection has to be carried out nu- 
merically. We calculated the projected mass density, the pro- 
jected velocity dispersion and a set of line profiles for our 
dusty galaxy models. In this section we compare the normal 
and the dusty projections of the kinematic quantities, and 
we investigate both the dependence upon the optical depth 
r and upon the dust geometry. 



4.1 Dependence upon optical depth 

In this section we keep the dust geometry fixed and vary the 
optical depth r. We adopt a modified Hubble profile {a — 3), 
which is slightly more extended than the stellar light. 



4-. 1.2 The projected dispersion 

The effect of dust on the projected dispersion profile cr^ {x) 
has a totally different nature than the effect on the projected 
density profile pp{x), since it is a normalized projection. For 
a particular line-of-sight x the weight function lC(x,r) is re- 
sponsible for the fact that not all parts of that line-of-sight 
contribute in the same proportion to the dusty projection, as 
is the case in the optically thin regime. From equation ( |l3[ ) 
and appendix U^ we can see that the weight function is a 
monotonically rising function of r for all lines-of-sight x, 
essentially because along a line-of-sight the near side dom- 
inates the light. Dust makes the contribution of the outer 
parts (in particular the near side) thus more important, rel- 
ative to the inner parts. 

In figure 0a,b we plot the projected velocity dispersion 
profile cTyix) for various optical depths. The explicit depen- 
dence of (J^{x) on r for a fixed line-of-sight x is shown in 
figure |3|c. From these figures it is clear that quantitatively, 
the effects of dust absorption on the projected velocity dis- 
persion are rather limited. Particularly for modest optical 
depths, r ^ 2, there is no major difference between normal 
and dusty projection of the dispersion. In figure Hd we plot 
the strength (in percentage) 



A<,(x) = 100 X 



o-pWl^^o" ol{x)\^ 



0-2(3;) 



(20) 
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Figure 2. The dependence of the projected density pp{x) upon the optical depth t. The dust model is the a = 3 model. In panel (a) 
the projected density is plotted as a function of the projected radius x for various values of r, which are indicated. In figure (b) we show 
the projected density explicitly as a function of the optical depth t, for three lines-of-sight x. The curve for a:: = can be calculated 
explicitly, as done in appendix Wi. The lower two panels show (c) the total integrated light Mobs, and (d) the core radius re and the 
effective radius re as a function of t. 



of the effect as a function of the anisotropy parameter q. It 
is shown for a set of lines-of-sight, with t — 2. The effect is 
only 1.6 per cent for the central line-of-sight, independent of 
the anisotropy of the model. At larger radii this anisotropy 
does become important, tangentially anisotropic systems are 
more sensitive to dust in their ap{x) profile than their radial 
counterparts. At x — Vc this drops to 1.2 per cent for the 
tangential and to 0.7 per cent for the radial model. 

Qualitatively, the general effect of dust is to lower the 
projected dispersion profile, and mostly Ao- will be a positive 
quantity. Only for the outer lines-of-sight of very radially 
anisotropic models the effect has the opposite sign, i.e. the 
projected dispersion is higher for the optically thick than 
for the optically thin case. In figure He we plot the critical 
line-of-sight Xc where Acr{xc) = as a function of q, for 
various values of r. Notice that its dependence upon the 
optical depth is very weak. 



4-1.3 The line profiles 

In figure ^ a number of normalized line-profiles lp(a;, v) are 
plotted, for a set of lines-of-sight varying from the centre to 
the outer parts of the galaxy. The line-of-sight velocities v 
are normalized with respect to the largest possible velocity 



one can observe at the particular line-of-sight x, being the 
escape velocity ^/2'ip{x), where ip{r) is the Plummer poten- 
tial (|7a|). The total area under the line profile, which equals 
the zeroth moment pp{x), is scaled to one. 

The effect of dust absorption on a normalized line- 
profile lp(a;, v) is at first order determined by the effect on its 
second moment, i.e. the projected dispersion a'^{x) at that 
projected radius. We can therefore copy most of our findings 
from the previous section. Only for higher optical depths 
and the inner lines-of-sight the infiuence of dust absorption 
is significant. In most cases dust absorption will turn the 
line profiles more peaked, and again tangential models are 
more sensitive to this than radial ones. 

A particular consequence is that the typical bimodal 
structure of the outer line profiles of tangential models can 
disappear (see the x — 1.5 line profile in figure Ml). Physically 
this can be understood in the following way. In a tangentially 
anisotropic model, most of the stars move on rather circular 
orbits, particularly in the outer parts. Since the density of 
the Plummer model follows a strong r~^ at large radii, the 
line profiles are primarily governed by the stars which move 
at orbits with radius r « a; (in both senses), causing the bi- 
modal structure. When the system is heavily obscured, the 
nearest parts of the line-of-sight contribute relatively more 
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Figure 3. The dependence of the projected dispersion o'p(x) upon the optical depth r. The dust model is the a = 3 model. In the 
upper two panels the projected dispersion is plotted as a function of the projected radius x for various values of r, which are indicated. 
It is shown (a) for a tangential model with q = —6 and (b) a radial model with q = 1. In panel (c) we show the projected dispersion 
explicitly as a function of the optical depth r, for four lines-of-sight x. It is plotted for the radial model q = 1. The curve for x = 
can be calculated explicitly, as done in appendix hI In figure (d) we plot the strength C\a(x) of the effect on the dispersion profile as a 
function of the anisotropy parameter q for three different lines-of-sight, which are indicated. The optical depth is t = 2. The last plot (e) 
displays the critical line-of-sight Xc for which ^aixc) = as a function of the anisotropy parameter q, for different values of r. Above 
these curves the parameter has negative values, below are positive values. 



than in the optically thin case. And these nearest parts con- 
tain stars moving on nearly circular orbits with larger radii, 
which have velocities that are gradually more perpendicular 
to the line-of-sight. Moreover these stars have smaller ve- 
locities anyway, since the circular velocity for the Plummer 
geometry declines as r~^" at large radii. Smaller line-of- 
sight velocities thus become predominant, which causes the 
bimodal structure to disappear. 



4.2 Dependence upon dust geometry 

In this section we vary the dust exponent between a = 1 
and a = 10, while the optical depth is held fixed. We choose 
large values for r (r = 10 and t = 20) so that the geometry 
effect is clearly visible. 



^.2.1 The projected density 

The dependence of Pp{x) on the dust exponent a is shown in 
figure 0a and pb, where we adopt the optical depth r — 10. 
As is to be expected, the projected density pp{x) is sensitive 



to variations of the dust geometry parameters. In the outer 
regions (the bottom curves in figure Hb) this is primarily due 
to the optical depth profile (114) which depends critically on 
a. Smaller values of a mean more dust and therefore more 
extinction. But also in the inner regions (black curve in fig- 
ure Hb) the extinction effect is larger for extended than for 
concentrated distributions, even though equal amounts of 
dust are present. At each projected radius x and for each 
optical depth r, the projected density Pp{x) is a rising func- 
tion of a, with lowermost value 

lim ppix) = i e-* ^ (21) 

a^i n (i + 2;2) 

As a consequence, the total integrated light Mobs increases 
with a too (figurepic), with exp (— §) for the overlying screen 
model being the minimum value. The more the dust is dis- 
tributed homogeneously over the line-of-sight, the more ef- 
ficient the extinction process. 

Also the morphology of the pp{x) profile is sensitive for 
variations of the dust geometry (figure Ha) . For small values 
of a the effect is small, and in the limit of the screen geom- 
etry the effect is nil : the projected density decreases with 
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Figure 4. A set of normalized line profiles for our standard dust model (ro = 1 and a = 3) for varying optical depth. The left, middle 
and right column correspond to respectively the radial, isotropic and tangential models. The line-of-sight x is indicated. Line profiles are 
shown for r = 0, t = 4, t = 10 and r = 20, going from black to light-coloured grey. Both the linc-of-sight velocities and the area under 
the curves are normalized, as described in the text. 



the same factor for all lines-of-sight. For larger values of a, 
the profile becomes flatter in tlie centre. For very condensed 
dust distributions the extinction is confined to the central 
parts of the galaxy. 

The dependence of the core radius re and the effective 
radius re on a is illustrated in figure pp. The maximum 
effect on these core size parameters is found at intermediate 
values of a. For very extended distributions on the one hand, 
the entire Pp{x) profile is affected with the same factor, and 
as a consequence the core size parameters are not affected. 
For very concentrated distributions on the other hand only 
the inner part of the core is affected, so that the total light 
output of the core remains nearly the same. 



4-2.2 The projected dispersion 

The primary effect of dust absorption on the projected dis- 
persion (Tp (x) at a certain line-of-sight x is to attach a weight 



to every part of that line-of-sight. Or, since both halves of 
our Plummer galaxy are dynamically equivalent, to every 
part of the nearer half of the line-of-sight. Since the param- 
eter Q controls the relative distribution of the dust along the 
lines-of-sight, this parameter determines in the first place the 
relative weights of the different parts. 

There are two extreme cases. For very concentrated 
models on the one hand, the entire first half of the line- 
of-sight is not affected. For very extended models on the 
other hand, the dust is effectively located between the ob- 
server and the galaxy, and all the parts of the line-of-sight 
X experience the same influence, and thus have the same 
relative weight. We obtain that for both cases, each part 
of the nearer half of the line-of-sight contributes equally to 
the projected dispersion, as in the optically thin limit. The 
dusty projected dispersion profile ap{x) therefore equals the 
optically thin profile (llq), no matter how large the optical 
depth r. 
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Figure 5. The dependence of the projected density Pp{x) upon the dust exponent a, with r = 10. In panel (a) the projected density is 
plotted as a function of the projected radius x for various values of a, which are indicated. Figure (b) shows the projected density as a 
function of a, for four hnes-of-sight x. The lower two panels show (c) the total integrated light Mobsi and (d) the core radius re and the 
effective radius re as a function of a. 



For the dust models considered here, the relative 
weights do differ from each other, and there is an effect on 
the projected dispersion profile. In the figures |6|a and |6p we 
plot the o'p{x) profile for various values of a, for a fixed op- 
tical depth T = 20. For every line-of-sight there is a critical 
exponent Oc that corresponds to the most efficient dust ge- 
ometry. This is illustrated in figure He, where a-p{x) is plotted 
as a function of a for various lines-of-sight. Notice that the 
position of the critical dust exponent Uc not only depends on 
the line-of-sight x, but also on the optical depth r and the 
anisotropy parameter q. This is shown in figure Hd, where a^ 
is plotted as a function of x for different values of r and q. 
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4-2.3 The line profiles 

The reasoning from the previous section does not only apply 
to the projected dispersion profile, but to all higher order 
moments of a line profile (for a definition, see appendix B). 
For each line-of-sight x, such a higher order projected mo- 
ment equals its optically thin value for very concentrated 
or very extended models, and there is an intermediate crit- 
ical value Qc for which the effect is maximal. This is illus- 
trated in figure M, where we plot the projected fourth mo- 



ment Tp{x) = iip{x)/pp{x) for various lines-of-sight, as a 



function of the dust exponent a. The analogy with figure 



dust exDonent a-1 

Figure 7. The projected fourth moment Tp(x) as a function of 
the dust exponent a for a set of lines-of-sight x. The optical depth 
is T = 20, the model is the radial q = 1 model. 



is obvious. Notice that the critical dust exponent Oc differs 
from moment to moment. For example, for the line-of-sight 
through the centre we find Qc — 2.115 for the projected 
dispersion, and Uc = 2.215 for the fourth moment. 

The same applies to the line profiles themselves (fig- 
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Figure 6. The dependence of the projected dispersion (Tp(x) upon the dust exponent a. The optical depth is fixed at the value t = 20 
for these plots. In the upper two panels o'p(a;) is plotted as a function of the projected radius x for various dust geometries, which are 
indicated. It is shown (a) for a tangential model with q = —6 and (b) a radial model with g = 1. In panel (c) we show the projected 
dispersion explicitly as a function of the optical depth a, for three lines-of-sight x. It is plotted for the radial model g = 1. In plot (d) we 
plot the critical dust exponent Oc (see text) as a function of the projected radius x for various optical depths and orbital models. The 
values of r are indicated. The solid lines correspond to the radial g = 1 model, the dotted lines to the tangential q = —6 models. 



ure pi). For the two limit cases, there is no difference be- 
tween optically thick and optically thin models, no matter 
how high the optical depth. And for each line-of-sight x, each 
optical depth r and each orbital model q there is a critical 
dust exponent Oc, for which the effect of dust absorption 
on the line profiles is maximal. The extent to which a dif- 
fers from this critical value determines the magnitude of the 
effect on the line profile. 



5 CONCLUSIONS 

Comparison of optical and IRAS data suggests that ellipti- 
cal galaxies are likely to contain a diffuse dust component, 
which is smoothly distributed over the galaxy (Goudfrooij 

6 de Jong 1995). Dust absorption affects the projection of 
all physical quantities along a line-of-sight. We constructed 
a four-parameter model to investigate the absorption effects 
on the projection of kinematic quantities. We compared the 
optically thick with the optically thin case, and investigated 
the effects of optical depth and dust geometry. 

We established that the projected mass density Pp{x) 
is substantially affected, even for modest optical depths 



(r ^ 2). Dust has also a non-negligible effect on the mor- 
phology of the pp{x) profile and on all global photometric 
quantities. On the other hand, the effects of dust absorption 
on the projected velocity dispersion profile ap{x) and on the 
line profiles are considerable only for higher optical depths 
(t <; 5). Moreover, the dust geometry is important ; both 
centrally concentrated and very extended dust distributions 
can remain invisible in the line profiles. 

Our models have some important applications, which 
are being investigated : 

• Due to the insensitivity of the Up (x) profile, dynamical 
mass estimates will hardly be affected by modest amounts 
of dust, while the light profile is sensitive for absorption. 
This should be taken in consideration in the interpretation 
of mass-to-light ratios in regions where no dark matter is 
expected : large mass-to-light ratios can indicate dust ab- 
sorption. 

• All global photometric properties of elliptical galaxies 
such as luminosities and radii are affected by dust. As WS96 
pointed out, this can lead to a systematic uncertainty in el- 
liptical galaxy distance estimation techniques that utilize 
these global properties, e.g. the L-a relation (Faber & Jack- 
son 1976) or the relations of the fundamental plane (Dressier 
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Figure 8. A set of normalized line profiles lp(a;, v) for various values of the dust exponent a and a fixed optical depth r = 20. The left, 
middle and right column correspond to respectively the radial, isotropic and tangential models. The line-of-sight x is indicated. Line 
profiles are shown for the optically thin model and for a = 5, a = 2 and a = 1.25, going from black to light-coloured grey. Both the 
line-of-sight velocities and the area under the curves are normalized. 



et al. 1987, Djorgovski & Davis 1987). With our models, we 
can not only quantify the effect of dust on the global pho- 
tometric properties, but also on the kinematics. 

• While large optical depths don't seem to be in accor- 
dance with global IRAS dust mass estimates, these models 
are nevertheless important on a local scale. Compact regions 
with high optical depths can be an explanation for the ob- 
served asymmetries and irregularities in kinematic data. In 
that case, one needs to combine (basically dust-free) near- 



infrared data with optical data. Such observations would 
help to constrain the diffuse dust masses in elliptical galax- 
ies, a problem which cannot be satisfyingly solved with the 
existing methods. 
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APPENDIX A: THE FUNCTION U^{z) 

For a > — 2 and for ^ 2: ^ 1 the function 11^(2:) is defined 



Uc.{z) 



2^2 



2 r 



'a+3) 



^/^^(^) J, 



(1-y')' dy. 

{l^y^f dy (Al) 



For natural values of a it can be easily calculated with the 
recursion formula 

\Jc{z) = U„-2(z) + -^ ^}J+2i Z {1-Z 



together with the seeds 

2 
U_i(z) — — arcsin(z) 

TV 

Uo(z) = z. 



2+1) 



(A2) 

(A3) 

(A4) 

The functions \Jaiz) can also be written in terms of the 
hypergeometric function, which is convenient for computa- 
tional purposes if a is non-integer, 

2 r(^) 



V^z) 



/^r(^; 



22 



fi(^ 



2 ' 2 ' 2 ' ^ y 



(A5) 



In the limit a —* —2 one finds the degenerate function 

,. J, , . r for < 2 < 1 
lim Uq(2) = < 
Q-.-2 ^ ' 1 for 2 = 1. 



(A6) 



APPENDIX B: THE MOMENTS OF THE 
CENTRAL LINE PROFILE 

The general moments ^i2n,2m,2iij') of a distribution function 
F{r,v) are defined as 



M2n,2m,2i(r) = j j j F [f , v) V^" Vg"^ vf dv , 



(Bl) 



where the integration ranges over all possible velocities. 
The moments fj,2n{r,x) of the projected velocity distribu- 
tion along a line-of-sight 2: at a particular point {x,r) in 
space are 



H2n{x,r) =^ . 



. , 2 2\^ — ^ 2i 

n\ [r - X ) X 



^J'2(n-i),2^,o{■r)■ (B2) 



To obtain the 2n'th moment of the line profile (or the pro- 
jected 2n'th moment) this has to be projected along the 
line-of-sight x, according to the recipe (|6a|,b), 



flp [X) 



IC{x,r] 



fi2n{x,r)rdr 



(B3) 



For the Plummer model, these projected moments can be 
evaluated analytically in the optically thin limit (Dejonghe 
1987), 



2n/\ 6-2 



r(6-g)r(n + i)r(^) 

r(n + 6-g)r(^) 



1 + x^ 



^ , q n+4 ^ n+5 x 

X3i^2(-n,-,^;l,^-,^^-^ 



(B4) 



yielding equation (hTI) and (hsl) for the cases n = and 
n = 1. 
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Table Bl. Analytical expressions for the central projected density Pp(0) and the central projected dispersion cr?^[0) as a function of the 
optic al de pth r. In all cases vq = 1, the value of a is indicated. The expressions for = 3 and a = 2 correspond to equations (B9|) 
and (plCJ). For a = 5 the projected density can be calculated by direct integration, for higher order moments no explicit expression can 
be deduced. 
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For dusty galaxy models, equation (B3) cannot be cal- 
culated analytically for general lines-of-sight.n Only for the 
central line-of-sight x = this is possible in a few cases, since 
the formulae become considerably easier. The sum (B2) re- 
duces to a single term. 



A*2n(0, r) = ^2n,o,o(r) 



3-2"-^ r(n-Ki) r(6~g) 
TTv^TT T {6 — q + n) 



(B5) 



(B6) 



yielding equation (iTq) and (fea) for n — and n — 1 respec- 
tively. If we substitute the weight function (hsl) of the King 
dust model, we obtain 



Mp"(0) = ce 
where we set 



cosh ( ^ Uo 



-3W)(1-. 



3.2"-i r(6-g)r(n-hi; 



tt-v/tt r (n -I- 6 — q) 



dz (B7) 



(B8) 



Since Uc,-3(2) can generally only be expressed as a hyperge- 
ometric function, we cannot evaluate this integral explicitly, 
except for some special values of the dust exponent a for 
which Uct-3(z) has a very simple form. For a = 3 the inte- 



gral in (B7) leads to a family of modified Bessel functions 
of the first kind (Gradshteyn & Ryzhik 1965, eq. 8.431), 

M^"(0)=4^cr(^) e-*r-^7^(^). (B9) 

For a = 2 we find (Gradshteyn & Ryzhik 1965, eq. 3.981) 

7r"+''cr(n + 4)e-5 cosh (f ) 



^p"(0) = <^ 



n^.:+''^'(r2 + (2i + l)27r2 
TT"+^cr{n + 4)e-i sinh(0 



^n 



(n + 3)/2 



+ (2j)27r2 



if n even 



if n odd. 



(BIO) 



From these expressions one can find the expressions for the 
projected density PpjO) and the projected dispersion crp(O), 
as tabulated in table Bl. 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science M^^X style file. 



* Except for the degenerate case a = 1, whe re /i2"(x) equals 
the prodiict of the optically thin moment (BJ) and the weight 
function (|l6|). 
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